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ABSTRACT 

We present a detailed study of the ionization fraction of the Barnard 68 
pre-stellar core, using millimeter H 13 CO + and DCO + lines observations. These 
observations are compared to the predictions of a radiative transfer model coupled 
to a chemical network that includes depletion on grains and gas phase deuterium 
fractionation. Together with previous observations and modelling of CO and 
isotopologues, our H 13 CO + and DCO + observations and modelling allow to place 
constraints on the metal abundance and the cosmic ionization rate. The H 13 CO + 
emission is well reproduced for metals abundances lower than 3 x and a 
standard cosmic ray ionization rate. However, the observations are also consistent 
with a complete depletion of metals, i.e. with cosmic rays as the only source of 
ionization at visual extinctions greater than a few A v . The DCO + emission is 
found to be dependent of the ortho to para H2 ratio, and indicates a ratio of 
~ 10~ 2 . The derived ionization fraction is about 5 x 10~ 9 with respect to H 
nuclei, which is about an order of magnitude lower than the one observed in the 
L1544 core. The corresponding ambipolar diffusion timescale is found to be an 
order of magnitude larger than the free fall timescale at the center of the core. 
The inferred metal abundance suggests that magnetically inactive regions (dead 
zones) are present in protostellar disks. 

Subject headings: astrochemistry - stars: formation — ISM: abundances — ISM: 
molecules — ISM: individual (Barnard 68) 



1. Introduction 



The ionization fraction (or the electron abundance) plays an important role in the chem- 
istry and the dynamics of prestellar cores. Because of the low temperature, the chemistry is 
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dominated by ion neutral reactions (IHerbst &: Klempererlll973l ). and electronic recombina- 
tion is one of the major destruction pathways for molecular io ns. Furthermore , the ionization 
fraction sets the coupling of the gas with the magnetic field (jShu et al.lll987l ). 



Several att empts have been made to estimate the electron fraction in dense clouds and 
prestellar cores jGuelin et alill982l : IWootten et allll982l ; Ide Boisanger et all 19961 ; I Williams et al 



1998; Caseiri et al. 



19981 ). These studies rely on measurements of the degree of deuterium 



fractionation (though the DCO + over HCO + abundance ratio for exam ple), which h as been 
found to be roughly inversely proportion al to the elec tron abundance (lLangerlll985l ). How- 
ever, this simple approach has caveats (jCasellil 120021 ) as it does not consider line-of-sight 
variations of the electron fraction. Large density gradients exist in prestellar cores, and 
therefore one may anticipate similar variations in the el ectron abundance. In addition, 
the f reeze-out of molecules onto the grain surfaces (e.g. iTafalla et al.l 120021 ; iBergin et al. 



2002h influence the degree of deuterium fractionation independently of the electron fraction 
( iCaselli et al.lll998l ). Finally, these studies usually consider simple chemical networks that 
may neglect important ingredients for the electron fraction. 

In this paper, we study the ionization fraction in the Barnard 68 core, using H 13 CO + 
and DCO + line observations. These observations are interpreted with a chemical network 
including gas-grain interactions that is coupled to a radiative transfer model. This technique 
allows us to infer the electron abundance along the line-of-sight, and to place constraints on 
the abundance of metals, the cosmic ray ionization rate, and the ionization state of material 
that is provided by infall to the forming proto-planetary disk. 

The paper is organized as follows: in Section [21 we present the observations. The model 
used to interpret these observations is detailed in Section [31 Implications of our findings are 
discussed in Section El and Section concludes this paper. 



2. Observations 

The H 13 CO+ (1-0) (u = 86.754288 GHz), and the DCO+ (2-1) (y = 144.077319 GHz) 
transitions were observed towards B68 (a = 17 h 22 m 38.2 s and 5 = -23°49'34.0"; J2000) in 
April 2002 and September 2002 using the IRAM-30m telescope. The core was mapped with 
a spatial sampling of 12". The half power beam size of the telescope is 29" at 87 GHz and 
17" at 144 GHz. System temperature were typically ~ 110 — 150 K at 3 mm and ~ 160 — 350 
K at 2 mm. Pointing was regularly checked using planets and was found to be better than 
~ 2". The data were calibrated in antenna temperature (T*) units using the chopper wheel 
method, and were converted to the main beam temperature scale (T mb ), using the telescope 
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efficiencies from the IRAM website. All observations were carried out in f requency switching 
mode. The H 13 CO + (1-0) data were also presented in lMaret et al.l (120061 ). 



Fig. [T]shows a comparison between the integrated lin e intensity maps o f DCO + (2-1) and 
H 13 CO+ ( 1-0) with the visual extinction map obtained bv klves et Zl (boOlh . The C 18 (1-0) 
map from lBergin et al.l (120021 ) is also shown. On this figure, we see that the peak of H 13 CO + 
(1-0) line emission does not correspond to the maximum visual extinction in the cor^j]. The 
C 18 (1-0) line emission shows a similar behavior: it peaks in a shell-like structure with a 
radius of ~ 50" around the maximum visual extinction. The DCO + (2-1) line emission, on 
the other hand, seems to correlate well with the visual extinction. These differences are likely 
a consequence of chemical effects. Because of the freeze-out on grain mantles, the abundance 
of CO and i ts isotopologues de crease by about two orders of magnitude towards the center 
of the core ( IBergin et al.ll2002l ). Since H 13 CO + is mainly formed from the reaction of 13 CO 



with H3 , its abundance is also expected to decrease towards the core center. DCO + should 
also be affected by the depletion of CO. However, the deuterium fractionation increases as 
CO is removed from the gas phase. Thus the disappearance of CO might be compensated by 
the increased deuterium fractionation. In the following, we interpret the emission of these 
species using a chemical model coupled with Monte-Carlo radiative transfer model, in order 
to derive precisely their abundance profiles. 



3. Analysis 



We have used a technique t hat combines the predict i ons of a chemical network with 
a Monte-Carlo radiative transfer (IBergin et al.ll2002l . 120061 ; iMaret et al.l 120061 ) . The outline 
of this technique is the following. Chemical abundances are computed as a function of the 
visual extinction in the core. Using these abundance profiles, the line emission is computed 
with a Monte Carlo radiative transfer code. The resulting map is convolved to the resolution 
of the telescope, and is compared to the observations. Free parameters of the chemical model 
(e.g. cosmic ionization rate, metal abundances, etc.) are adjusted until a good agreement is 
obtained between the model and the observations. Thus, this technique allows for a direct 
comparison between the predictions of the chemical network and the observations. 



We have used the chemical network of IBergin et al.l (119951 ). This network is contains 
about 150 species (including isotopologues, see below), and focuses on the formation of 



1 This is also clearly seen on Fig. [21 which shows the H 13 CO + (1-0) line emission as a function of the 
visual extinction. The line emission increases as a function of the A v between and 20, but decreases at 
A v ~ 20. 
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simple molecules and ions (e.g. CO and HCO + ). The network includes the effect of deple- 
tion on grains, and the desorption by thermal evap oration, UV photons, and cosmic rays 
( IHasegawa fc Herbstl Il993l ; iBringa &: Johnson! 120041 ) . It also includes the e ffect o f fraction- 
ation of 13 C and 18 0, using the formalism described by lLanger &: Penziad (119931 ) . We have 
extende d this network to inc lude the effect of deuterium fractionation, following the approach 
used by iMillar et al.l (Il989l ). Because of the importance of multiply deuterated species in 
the de uterium fractionation process, these species were also included in the network, fol- 
lowing |Roj2exts_el_aL| (120041 ) . It also include neutralization reactions of ions on negatively 
charged grains. T he predictions of our network were checked against the UMIST network 
( IMillar et al.lll997l ) for consistency. 



We adopt the density profile determined by lAlves et al.l (120011 ). from observations of 
near infrared extinction from background stars. This profile is assumed to be constant as a 
funct ion of time. The dus t temperature profile was computed using the analytical formulae 
fro mlZucconi et al.l (120011) . For the gas temperature we have adopted the profile determined 
bv iBergin et all (bood l from observations and modelling of CO and its isotopologues. The 
gas temperature is relatively low (7-8 K), and increa ses sightly (10- 11 K) at the center of 
the core as indicated by ammonia lines observations (ILai et al.l 120031 ). This increase in the 
temperature is a result of grain coagulation at the center of the core, which produces a 
thermal decoupling between the gas and the cooler dust. 

The cloud is supposed to have the initial composition summarized in Table [TJ In our 
model, we assume that the density profile of the core does not evolve with time. Therefore, 
we also assume that the chemistry has already evolved to a point where hydrogen is fully 
molecular, and all the carbon i s locked into CO. Our treatment of the initial atomic oxygen 
pool deserves special mention. IBergin fc Snelll (120021 ) examined this question in the context 
of the non-detection of water vapor emission in B68 by SWAS. They found that if atomic 
oxygen were present in the gas phase in the dense core center, then the well studied reaction 
chain that forms H2O (via H30 + ) would have yielded detectable water vapor emission. The 
simplest way to stop this reaction chain is to remove the fuel for the gas-phase chemistry: 
atomic ox ygen. This happen s when oxygen is trapped on grain surfaces in the form of water 
ice (e.g. IBergin et al.ll2000l ). Thus we have assumed initial conditions in which all non- 



refractory oxygen is in the form of water ice and CO gas with no atomic oxygen left. In this 
fashion our initial abundances assume the core formed out of gas that reached at A v ~ 2 
- where H2 and CO have formed and water ice mantles are observe d. On the other hand, 
nitrogen is assumed to be mostly in atomic form (IMaret et al.l 120061 ) . 



A grain size of 0.1 fim is assumed. The cosmic ray ionization rate and the abundance 
of low ionization potentials metals (< 13.6 eV) are free parameters of our study (see §4.11 
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and §4.2p . In our models we combine all metals (e.g. Fe + , Mg + , ...) into one species, labeled 
as M+ with the Fe+ recombination rate of a(M+) = 3.7 x 10~ 12 (T/300 K)~ a65 cm 3 s -1 . Due 
the low ionization potential these metals are assumed to be fully ionized at the start of the 
calculation. The network also includes the neutralization of ions of negatively charged grains 
with one electron per grain. 



The core is assumed to be b athed in a UV field of 0.2 (in Habing units; Il968l ). as 
determined by iBergin et al.l (120061 ) . The chemical abundances are co mputed as a functio n 
of time by solving the rate equations using the DVODE algorithm (IBrown et al.l 119890 . 
This is done until a time of 10 5 yr is reached. This corresp onds to the "best-fit" model of 
Bergin et al.l ( 20061 ). However, as discussed by IBergin et al.l . this time is a lower limit of the 
real age of the cloud, since the CO is assumed to be pre-existing at t — in these models. 

Modeling the line emission requires the knowledge of yelocity profile in the core. As a 
first approach, we have neglected systematic motions (see lLada et al.l 120031; iRedman et al. 
2006 ). and we have used the turbulent velocity profile determined by Bergin et al. ( 20061 ) 



from C 18 and 13 CO lines. The turbulent velocity is ~ 0.3 km/s at the edge of the cloud, 
and decreases significantly (~ 0.15 km/s) towards the center of the core. 



4. Results 



4.1. Metals depletion 



Metals ions (e.g. Fe + and Mg + ) play an important role in in setting the electron abun- 
dance in pre-stellar cores, because they are destroyed relatively slowly by radiative recombi- 
nation. For example, the recombination rate of is four order of magnitude higher than 
the rate for Fe + . 



Guelin et al.l (119821 ) measured the electron abundance in a sample of dense molecular 
clouds using HCO + and DCO + line observations, and obtained values comprised between 
10~ 8 and 10~ 7 . The authors concluded that the metal abundance is lower than 10 -7 in 



these clouds. ICaselli et al.l (119981 ) determined the electron abundance in a sample of twenty 
four low-mass isolated cores (with embedded stars and starless - similar in properties to 
B68) from CO, HC O + and DCO + observations, and obtained values in the range 10~ 8 - 



10~ 6 . ICaselli et al.l argued that the differences between cores are due to changes in metal 
abundance and a variable cosmic ionization rate ((). The best fit between their chemical 
model pr edictions and the obs ervations indicates metal abundances in the range 2 x 10~ 9 — 
3 x 10~ 7 . IWilliams et al.l (119981 ) determined the electron abundance in a similar sample of low 



mass cores using a slightly different approach, and obtained metal abundances comprised 
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between 5 x 10~ 9 and 4 x 10~ 8 (assuming a constant (). All these studies indicate low 
metal abundances with respect to their solar values.. Indeed, observations of FUV Fell 
absorption lines, and other metal lines, towards diffus e clouds find depletion factors of over 



two order of magnitud e with respect to solar values (jSavage &: Bohlinl Il979l ; IJenkins et al. 
19861 : ISnow et alJbooj l. 



Our H 13 CO + (1 — 0) observations can be used to set limits on the metal ion abundance 
in B68. H 13 CO + is sensitive to the electron abundance inside the core, because it is mainly 
destroyed by electronic recombination. It is also sensitive to the H3 and 13 CO abundances, 
since it is formed from the reaction between these two species. H3 itself is mainly formed 
from H2 ionization by cosmic rays. The remaining parameter in determining the chemi- 
cal abundance profil e is the time depen dence of the chemistry. In this case, our analysis 
is simplified because iBergin et al.l (120061 ) used multiple transitions of 13 CO and C 18 and 
a similar modeling technique to derive the 13 CO abundance and constrain the "chemical 



ag 
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of Barnard 68 to t = 10 5 yr. Thus, the only free parameters for our modeling of the 
H 13 CO + emission are the cosmic ionization rate £ and t he metal ion abund ance. These two 
parameters are difficult to constrain simultaneously. In iMaret et al.l (120061 ). we found that 
the H 13 CO + (1 — 0) line emission in B68 is well reproduced by our chemical network if one 
assume a metal abundance of 3 x 10~ 9 with respect to H nuclei and a standard cosmic ion- 
ization rate (£ = 3 x 10 -17 s -1 , see next section). In the following, we explore the parameter 
space into more details to place constrains on the metal abundance in the core. 

On Fig. [21 we show the predicted intensity of the H 13 CO + (1-0) line for different metal 
ion abundances and cosmic ionization rates. In these models, metals are assumed to be 
initially fully ionized. In Fig. [2], we see that for ( = 3 x 10~ 17 s _1 , our model predicts the 
same intensities for x(M + ) = and x(M + ) = 3 x 10~ 10 . The predicted emission is in fairly 
good agreement with the observations. On the other hand, for a higher metal abundance 
(x(M + ) = 3 x 10~ 9 ) the model predicts a intensity slightly lower than the observed, but is 
in better agreement with the observations at the center of the core. A metal abundance of 
3 x 10~ 8 is clearly ruled out by the model and observation compariso n. We conclude that 
x(M + ) < 3 x 10~ 9 . This v alue is at the low end of the one obtained by lCaselli et al.l (119981 ) 



and I Williams et al.l (119981). Compared to the abundance of metals in the solar photosphere 
[x(M) ~ 8.5 x 10~ 5 ; I Anders &: Grevesselll989l ). this represent a depletion factor of more than 



four orders of magnitude. Indeed, our observations are also fully consistent with a complete 
depletion of metals in the core, i.e. with cosmic rays as the only source of ionization at A v 
greater than a few magnitudes (see Fig. [2]). It should be noted, however, that result depends 



See discussion in IBergin et al.l (|2006T ) on the meaning of this "chemical age" . 
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on the value of ( adopted. For example, our observations are fully consistent with a cosmic 
ionization rate of 3 x 10~ 16 s _1 and x(M + ) = 3 x 10~ 8 . The effects of varying ( are discussed 
in the next section. 



4.2. Cosmic Ray Ionization Rate 



Cosmic rays play a crucial role in the chemistry of pre-stellar cores, because they set the 
abundance of the pivotal Hjj" ion, and are the only source of ionization at A v greater than a few 
mag nitudes. Despite of it s importance, the cosmic ray i onization rate is difficult to constrain 



see 



Le Petit et al.l 12004 . Ivan der Tak et al.l 120061 and iDalgarnol 120061 for recent reviews 



7 x lO^s" 1 



Early estimates in diffuse clo uds from HD and OD observations indicate ( 
( Ivan Dishoeck fc Black! Il986l ). a value in agree ment with the lower limit of 3 x 10 _17 s _1 
measured by the Voyager and Pioneer satellites (jWebberl 119981 ). Hg ob servations towards C 



Persei cl oud suggest a s i gnific antly higher rate = 1.2 x 10~ 15 s -1 ; iMcCall et al.ll2003l ). 
However, iLe Petit et al.l (120041 ) argued that a value of ( = 2.5 x 10 -16 s _1 is more consistent 
with both Hjj" and HD observations. In de nser regions, HCO + observat i ons i ndicates a 



lower ionization rate than in diffuse clouds: Ivan der Tak fc van Dishoeck! (120001 ) obtained 
( = (2.6 ± 1. 8) x 10~ 17 s~ 1 f rom H CO + line observations towards massive protostars. In pre- 
stellar cores, ICaselli et al.l (119981 ) inferred a value comprised between 10 -18 and 10 -16 s _1 . 
The difference in the cosmic ray io nization rate between d iffuse and dense clouds could be 
due to the scattering of cosmic rays (IPadoan fc Scaloll2005r). In addition, l arge variations are 



inferr ed as a function of the Galactic Center distance (jOka et al.l 120051 ; Ivan der Tak et al. 
20D6h . 



Cosmic rays are also heating agents of the gas. iBergin et al.l (120061 ) examined the value 
of ( in B68 by comparing th e predictions of a chemical and thermal model to observations 
of CO and its isotopologues. IBergin et al.l found that their model provide reasonable fits to 
the data for ( = 1 - 6 x 10~ 17 s _1 . Their "best fit" model has C = 1.5 - 3 x 10~ 17 s -1 . Here 
we examine the constraints placed by our H 13 CO + observations. On Fig. [2J we see that our 
model produces a good fit to the data for ( = 3 x 10~ 17 s^ 1 , except for x(M + ) = 3 x 10~ 8 , 
where the model predictions underestimate the observation by a factor two. Models with 
( = 3 x 10~ 18 s _1 , consistently underestimate the observations. Conversely models with 
( = 3x 10~ 16 s " 1 overestima t e the model, except the one with x(M + ) < 3 x 10~ 8 . This in 
agreement with IBergin et al.l (120061 ). who found that their observations are not reproduced 
by models with ( = 6x 10 -16 s _1 . 



To summarize our conclusions regarding the metals abundances and the cosmic ioniza- 
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tion rate, models with x(M + ) < 3 x 10 9 provide a good agreement with the data§, although 
the model with x(M+) = 3 x 1(T 8 and ( 



3 x 10 16 s 1 is also consistent with our data. 



However values of ( greater that 6 x 10 _17 s _1 are ruled out by lBergin et al.l (120061 ) based on 
core thermal balance. On the other hand, models with ( = 3 x 10~ 18 s _1 always underesti- 
mate our observations. We conclude that ( = 1 — 6 x 10~ 17 s -1 , and x(M + ) < 3 x 10~ 9 in B68. 
This implies that the abundance of ionized metals is reduced in the center of B68. Charge 
transfer from molecular ions (e.g. H^, HCO + ) to metals can be important and a reduction 
in the abundance of ionized metals also re quires lowering the neutral metal abundance. In 
the case of Fe a p otential reservoir is FeS (Keller et al.ll2002l ). or organometallic molecules 
( ISerra et al.lll992l ). An other possility is that Fe is incorporated into grain cores. 



4.3. Ortho to para H 2 ratio 

The ortho to para H 2 ratio influences the degree of ion and molecule deuteration in 
prestellar cores (IPineau des Forets et al.l Il99ll ; iFlower et al.l l2006bl). In the gas phase, deu- 



Roberts et al. 2004, and 



terium fractionation is mainly due to the following reaction (see 
references therein): 

H+ + HD ^ H 2 D+ + H 2 (1) 
The reverse reaction has an activation barrier of ~ 232 K and therefore the reaction becomes 



essentially irreversible at low temperature. iGerlich et al.l (120021 ) measured the forward and 
reverse rates of the above reaction at 10 K, and found them to be very different than com- 
monly adopted values. The forward reacti on rate was found to be about five times higher 
than previous estimates (jSidhu et al.lll992l ). while the reve rse reaction rate wa s found to be 



five orders of magn itude larger than previously used (e.g. iCaselh' et al.lll998l ). In addition, 
Gerlich et al.l (120021 ) determined via a laboratory measurement that the reverse reaction rate 
is very sensitive to the ratio of ortho to para molecular hydrogen. This is because o-H 2 , in 
its ground rotational level (J = 1) has an higher energy (AE ~ 170.5 K) when compared 
to the ground state of p-H 2 (J = 0). Consequently, o-H 2 can more easily cross the energy 
barrier than p-H 2 , and the rate of the reverse reaction increases with the ortho to para H 2 
ratio. 

Our DCO + observations can be used to estimate the H 2 D + abundance, and thus the 
efficiency of the deuterium fractionation process. DCO + is mainly formed by the following 



3 For simplicity, we have assumed that the initial x(M + ) abundance is constant as a function of a radius. 
A better fit to the observations might be obtained with a variation of x(M + ) with the radius. 
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reaction: 

H 2 D+ + CO -> DCO+ + H 2 (2) 

and is mainly destroyed by electronic recombination. Thus the DCO + emission depends on 
both the CO and H2D + abundances, the electron fraction (induced by cosmic-rays and by 
pre-existing metal ions), the ortho to para H2 ratio, and on time. Here we benefit from our 
previous analysis of CO which constrained the CO abundance and "chemical age" and our 
analysis of H 13 CO + which limit the metal ion abundance and cosmic ray ionization rate. 
Thus the primary free parameter is the ortho to para H2 ratio (o/p) when we adopt our best 
fit parameters of x(M + ) = 3 x 1(T 9 and ( = 3 x 10 -17 s -1 . 

On Fig. [3j we compare the observed DCO + (1-0) line emission as a function of A v , with 
the predictions of our model for different ortho to para H 2 ratio. Note that in these models, 
no o/p conversion is considered: the ortho to para H 2 ratio is assumed to be constant. The 
best agreement^ between the observations and the model is obtained for an o/p ratio of 
~ 1.5 x 10 -2 , well above the Boltzmann equilibrium value at 10 K (3.5 x 10~ 7 ). Fig. [3] also 
show the derived DCO + abundance inside the core. The abundance peaks at an A v of ~ 5, 
and decreases slightly towards the core center, as a consequence of CO depletion (see Section 

EI]). 

It is interesting to c o mpar e the o/p H 2 ratio we obtain with the predictions of other 



models. I Walmsley et al.l (120041 ) modeled the o/p H 2 ratio in prestellar cores, assuming a 
complete depletion of heavy elements. In their model, an initial o/p H 2 ratio of 3.5 x 10~ 7 
is assumed. For a density of 10 6 cm -3 , steady-state is reached in 10 5 yr, a time comparabl e 



to the age of B68 inferred from CO depletion observation and modeling flBergin et al.ll2006l ). 



At steady state, the o/p H 2 ratio obtained is 6 x 10 5 , i.e. about tw o orders of m a gnitud e 



lower than the value determined in this work. However, as noted by iFlower et al.l (l2006al ). 
the ortho to para H 2 ratio conversion reactions are very slow, and it is not clear if the steady 
state equilibrium is reached in molecular clouds prior to the formation of dense c ores. Using a 



i nitial o rtho to para ratio of 3 (a value appropriate for H 2 formation on grains) . IFlower et al. 



(j2006bl ) obtain a steady state ratio of 3 x 10~ 3 . This value, although still about a factor 5 
lower, is in better agreement with our estimate. We note that for o/p = 3 x 10 -3 , our model 
predicts an DCO + (1-0) emission about 2 times higher than the observations (see Fig. [3D- 



4 Although the model predicts the correct intensity at the core center, one can note that emission at lower 
A v is slightly underestimated. This may suggest an H2 o/p variation with the radius of the core: increasing 
the ratio at low A v would increase the DCO + abundance emission in this region and would probably produce 
a better fit. 
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Discussion 



5.1. Electron abundance and main charge carriers 



On Fig. m we show the derived electron and main ions abundances inside the core. 
The electron abundance is ~ 5 x 10~ 9 with respect to H nuclei throughout most part of 
the core. At low A v , the electron abundance increases as a result of photo-dissociation of 
CO. In this region, the most abundant ion is C + . At higher A v , the most abundant ion 
is H^, which caries about ~ 20% of the electric charge. The remainder of the charge is 
shared between more complex ions. Deuterated ions do not contribute significantly to the 
ionization fraction. In the innermost region of the core, where the deuteration increases as 
a result of CO depletion, the main deuterated ion, Djj~, is about ten times less abundant 
than Hjj\ H 2 D + and D 2 H + have simila r abundances (2 x 10~ n with respect to H). This is 



in agreement with recent observations (jVastel et al.1 12004 ) . 



Recently, iHogerheijde et al.l (120061 ) reported a probable detection of o — H 2 D + funda- 
mental line towards B68 which can be compared to our model predictions. The measured flux 
is however quite uncertain, given the relatively low signal to noise ratio of this observation 
(2.7a and 5.2a on the peak and integrated i ntensity, respective ly). Assuming a thermal exci- 
tation (10 K) and optically thin conditions, IHogerheijde et al.l derive a H2D" 1 " column de nsity 
of 1.5 x 10 12 cm -2 . Assuming a H2 column density of 3.6 x 10 22 cm~ 2 ( lAlves et al.ll200ll ). this 
corresponds to an H 2 D + abundance of 2.1 x 10~ n with respect to H nuclei, averaged in the 
APEX beam (17"), with respect to H nuclei. This is in excellent agreement with our model, 
which predicts an H 2 D + abundance of 2 x 10~ n , roughly constant across the envelope. Of 
course, if the excitation is non-the rmal, the detection implies an higher abundance. Assum- 
ing a 5 K excitation temperature, IHogerheijde et al.l derive a beam averaged abundance of 
1.5 x 10~ 10 with respect to H nuclei. This is about an order of magnitude higher than our 
model predictions. Since no collisional rates exist in the literature for H 2 D + , it is unclear 
whether or not the excitation of this line is thermal. IHogerheijde et al.l estimate a critical 
density of 2 x 10 6 cm -3 , which exceeds the density at the center of B68 (3 x 10 5 cm -3 ) 
by about an order of magnitude. However, the collisional rate, and therefore the critical 



densi ty, is uncertain by an order of magnitude (Ivan der Tak et al.1 120051 ; IHogerheijde et al. 



20061 ). Our model predictions regarding the deuterium chemistry could be also tested via 
observations of the D 2 H + li o — lo,i [y = 691.66044 GHz). Assuming a excitation temper- 
ature of 10 K, we predict a line intensity of 10 mK. Unfortunately, this is too weak to be 
detected with current ground based telescopes. 

We w ould like to c o mpar e the electron abu ndance profil e we obtained with the one 
derived by ICaselli et al.1 (120021 ) in L1544. In the ICaselli et al.l best fit model, the electron 
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abundance at the center of L1544 is 5 x 10~ 10 (with respect to H), while we obtain an 
electron abundance an order of magnitude higher at the center of B68. These differences 
are probably a consequence of different central densities: the L1544 central density is about 
an order of magnitude higher than the one of B68, and the electron fraction is expected 
to scale as nfH^) - 1 / 2 ( McKed Il989f ) . Another important difference is the dominant ion: 
Caselli et al.l (120021 ) predicts that the most abundant ion is H 3 + , while in our modeling 



main charge carrier is Hj~. These differe nces are due to different assumptions on the atomic 
oxygen abundance. ICaselli et al.l (120021 ) assumes that oxygen is initially mostly atomic. As 
a consequence, the H 3 + abundance is relatively large, because atomic oxygen reacts with 
to form H 3 + (after successive protonations by H 2 followed by recombination). In our 
modeling, oxygen is assumed to be initially locked in water ices and gas phase CO (see Table 
[1]), and the atomic oxygen gas phase abundance is relatively low. 

Finally, we would like to commen t on th e effect of grain size evolution on the electron 
fraction in the core. IWalmsley et al.l (120041 ) computed the electron abundance and main 
charge carrier in a prestellar core for different grain sizes. For a grain size of 0.02 /j,m, 
the main charge carrier in their model is H^, while for larger grains (0.1 fim), the most 
abundant ion becomes H + . In their models, H + recombines primarily on grains, while 
recombines with free electrons. Since the recombination timescale on grains depends on the 
grain size, the H + over abundance ratio, and in turn the electron abundance, depends 
on the grain size as well. However, these models assume a complete depletion of heavy 
elements, which is not the case for B68. In B68 we do find evidence for strong molecular, 
but not complete, heavy element freeze-out, at the core center. The reaction with H + with 
molecules containing these elements (e.g. NH3, OH, ...) can transfer the charge to molecular 
ions with faster recombination timescales. This would probably reduce the dependence of 
the electron abundance on the grain size. 



5.2. Core stability 



The electron abundance in the core is also important for its dynamical evolution, since 
its affects the efficiency of ambipolar diffusion. In a weakly ionized sub-critical core, the 
ions are supported again st collapse by th e magnetic field, but neutrals can slowly drift with 
respect to the ions fsee IShu et al.l 119871 . for a review). The timescale for this phenomenon 
is given by IWalmsley et al.l (120041 ) : 



Tad 



TiGm 2 n 



11; 



m~m„ 



n n rrii + m. 



■{crv)i 



(3) 
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where G is the gravitational constant, m n and im are the masses of the neutrals and the 
ions respectively, n n and n^ are the number densities, (<Jv)i n is the rate coefficient for the 
momentum transfer, and the summa tion goes ove r all ions. At low temperature, the rate 
coefficient for momentum transfer is ([Flower! l2000l ) : 



(<rv) ir . 



2ne I a 



rrii + m r 



rrumr, 



1/2 



where a is the polarizability of H2. Assuming that Hjj" is the dominant ion, we obtain: 



r ad ~2xl0 14 x(e) yr 



(4) 



(5) 



where x(e) is the electron abundance, with respect to H. Thus at the center at the core, the 
ambipolar diffusion timescale is 10 6 yr. It is interesting to compare this to the free fall time 
scale, which is given by: 

where p = n H2 mH 2 is the mass density. When expressed as a function of n H2 , this gives: 



r s = 3.6 x 10 7 n H l /2 yr 



(7) 



At the center of B68 we obtain rg = 7 x 10 4 yr which is about an order of magnitude faster 
than the ambipolar diffusion timescale. Thus, if present, the magnetic field may provide an 
important source of support. 

The strength of th e magnetic field that is neede d to support the cloud can be obtained 
from the critical mass (jMouschovias & Spitzerl 119761 ): 



M 



0.13 



where <$b = ttR 2 B is the magnetic flux, R is the core radius, and B is the magnetic field 
strength. The strength of the magnetic field that is needed to support the cloud is therefore: 



B 



G l ' 2 M 
" 0.13 tiR 2 

12500 AU and M 



(9) 



where M is the mass the core. Using R = 12500 AU and M = 2.1 M Q flAlves et ajj|200l| ), we 
obtain a critical magnetic field of 76 /iG for B68. No magnetic field measurements for B68 
exist in the literature, but we can compare this value to th e one m easu red in other cor e s from 
dust sub-millimeter polarization. I Ward- Thompson et al.l (120001 ) and ICrutcher et al.l (120041 ) 
measured plane- of-the-sky magne tic field strengths of 80 /xG in L183, 140 /iG in L1544 and 
160 fiG in L43. iKirk et al.l (120061 ) measured lower fields of 10 and 30 /iG in the L1498 and 



13 



L1517B cores. Therefore, if the magnetic field strength in B68 is at the lower end of the 
values measured in other cores, then it might be super-critical (i.e. the magnetic field is too 
weak to balance gravity). If it is higher, then the core is probably sub-critical. One may 
argue B68 has nearly round shape (albeit with an asymmetrical extension to the southeast), 
which potentially is indicative of a weak magnetic field. 



5.3. Implications of the metals depletion for accretion in protostellar disks 



One important conclusion of this study is the large metal depletion inferred for B68. 
Here we examine the implication of this findings for the mechanism of angular momentum 
transport in protostellar disks. The most favored theory for angular momen tum transport in 
disks predicts that accretion occurs via magneto-rotational instability (MRI; iBalbus fc Hawley 
19911 ) which produces MHD turbulence. Since this is a magnetic process, the ion-neutral 
coupling is therefore important. Typically, the ionization f raction should be gre ater than 



10 12 for disks t o be able to sus tain MHD turbulence (see lllgner fc Nelson! 12006k and ref- 



erences therein). iGammid ( 119961 ) suggested a model wherein the accretion is layered. The 
electron abundance is high at the surface of the disk, because of the ionization of the gas 
by UV, X-rays and cosmic rays, but it decreases towards the mid-plane. Thus disks may 
have a magnetically active zones at high altitude, where the electron fraction is sufficient 
to maintain MHD turbulence, and dead zones, closer the mid plane of the disk, where the 
electron fraction is lower, and accretion cannot occur. 

Our results have some import on this process because the chemical structure of the 
pre-stellar stage sets the initial chemical conditions of the gas that feeds the forming proto- 
planetary disk. Because of their influence on the ionization fraction, metal ions can have 
dramatic eff ects on the size of the dead zone, assumin g that they are provided by infall 
to the disk (IFromang et al.l 120021 ; lllgner fc Nelson! 120061 ) . The latter authors computed the 
ionization fraction in a protostellar disk, and found that for x(M + ) < 3 x 10 -10 , the dead 
zone extend between 0.5 and 2 AU, while it disappears completely for x(M + ) > 10 -8 . In 
B68, we obtain a metal abundance of x(M + ) < 3 x 10~ 9 , which is below the threshold for 
a complete disappearance of the dead zone. Thus, if B68 is representative of the initial 
conditions for the formation of protostellar disks, and cosmic rays do not penetrate deeply 
to the midplane, dead zones should exist in those disks. 
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6. Conclusions 

We have presented a detailed analysis of the electron abundance in the B68 prestellar 
core using H 13 CO + (l — 0) and DCO + (2 — 1) line observations. These observations were 
compared to the predictions of time dependent chemical model coupled with a Monte-Carlo 
radiative transfer code. This technique allows for a direct comparison between chemical 
model predictions and observed line intensities as a function of radius (or the visual extinc- 
tion) of the core. Our main conclusions are: 

1. The metal abundance is difficult to constrain independently from the cosmic ioniza- 
tion rate. However, accounting for thermal balance considerations and to reproduce 
H 13 CO+(l - 0) emission we estimate that x(M+) < 3 x KT 9 and C = 1 - 6 x KT 17 s -1 . 

2. The DCO + (2 — 1) line emission is sensitive to the ortho to para ratio. The emission is 
well reproduced by our model for an ortho to para ratio of 1.5 x 10~ 2 , well below the 
equilibrium value, and in reasonable agreement with previous work. 

3. The inferred electron abundance is 5x 1CT 9 (with respect to H), and is roughly constant 
in the core at A v > 5. It increases at lower A v because of the photo-dissociation of 
CO and photo-ionization of C. In the dense part of the core, the dominant ion is 
H^. H 2 D + and D 2 H + have similar abundances and are about two of magnitude less 
abundant than Hj. In the center of the core, our model predicts Djj~ to be the most 
abundant deuterated ion. 

4. The inferred electron abundance implies an ambipolar diffusion timescale of 10 6 yr at 
the center of the core, which is about an order of magnitude higher than the free fall 
timescale (rj = 7 x 10 4 yr). 

5. The metal abundance we obtain is below the threshold for protostellar disk to be fully 
active. Consequently, if the chemical composition of B68 is reprentative of the initial 
conditions for the formation of a disk and cosmic rays do not penetrate to the disk 
mid-plane, then dead zones should exist in protostellar disks. 

Both authors are grateful to C. Lada for a fruitful collaboration that led to this work 
and to T. Huard and E. Aguti for obtaining a portion of these data. We are also grateful to 
the referee and to the editor J. Black for useful and constructive comments. S. M. wishes to 
thanks H. Roberts for helping us testing the predictions of our network for deuterated species, 
E. Herbst for useful discussions about chemical reactions rates, E. Roueff and M. Walmsley 
for discussions about the ortho to para ratio, L. Hartmann and F. Heitch for discussions 
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Fig. 1. — Compari son between integrated intensity maps (contours) of C 18 (1-0) (left, from 
Bergin et al.ll2002l ). H 13 C0 + (1- 0) (center) and DC O + (2-1) (right) superposed on the map 



of visual extinction obtained by lAlves et all (120011 ). C 18 (1-0) contours start at 0.2 K km 
s _1 and step by 0.2 K km s _1 . H 13 (1-0) contours start at 0.15 K km s _1 and step by 0.15 
K km s . DCO + (1-0) contours start at 0.1 K km s _1 and step by 0.1 K km s -1 . The Ay 
image range from to 27 mag. 
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Fig. 2. — Comparison between the observations and the model predictions of the H 13 CO + 
(1-0) line for different cosmic ionization rates (£) and metal abundances. Blue points with 
error bars (1 a) represent the observed integrated line intensity as a function of the visual 
extinction in the core (A v ). Dashed, solid and dash-dotted lines represent the model predic- 
tions for ( = 3 x 10~ 18 , 3 x 10 -17 and 3 x 10 _16 s _1 respectively. In panel (a) a complete 
depletion of metals is assumed. In panels (b), (c), and (d), a metals abundance of respectively 
3 x 10- 10 , 3 x 10~ 9 and 3 x 10" 8 is assumed. 
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Fig. 3. — Comparison between the observations and the model predictions of the DCO + 
(2-1) line for different ortho to para H 2 ratios. Blue points with error bars (1 a) represent 
the observed integrated line intensity as a function of the visual extinction in the core (A v ). 
The dashed, dotted and solid black lines show the predicted line intensity for different ortho 
to para H 2 ratios. The inset show the corresponding abundances of DCO + . 
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Fig. 4. — Derived abundances inside B68 for the electrons and main ions. Abundances are 
relative to H nuclei. 
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Table I. Initial abundances. 



Species Abundance 3, 



H 2 


0.5 




He 


0.14 




H20i C es 


2.2 x 10~ 


-4 


- L1 2 w ices 


4.4 x 10" 


-7 


CO 


8.5 x 10" 


-5 


13 co 


9.5 x 10- 


-7 


c 18 o 


1.7 x 10" 


-7 


N 


1.50 x 10 


-5 


N 2 


2.5 x 10- 


-6 


HD 


1.6 x 10" 


-5 


Grains 


io- 12 





a Relative to H nuclei. 



